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Abstract 

We find the time evolution of the system of two non-interacting unstable particles, 

■ distinguishable as well as identical ones, in arbitrary reference frame having only 
\^ . the Kraus operators governing the evolution of its components in the rest frame. We 

than calculate in the rigorous way Einstein-Podolsky-Rosen quantum correlation 

■ functions for K^K^ system in the singlet state taking into account CP- violation 
\^ . and decoherence and show that the results are exactly the same despite the fact we 

I treat kaons as distinguishable or identical particles which means that the statistics 

of the particles plays no role, at least in considered cases. 

Dh: 

1 Introduction 

. 5t , In the recent years the possibility of testing Bell-CHSH inequalities [1,2] in 

the system of correlated neutral mesons has attracted some attention (see e.g. 
[3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19]), because kaons and 5-mesons are 
detected with higher efficiency than photons which maybe allows one to close 
the detection loophole (see e.g. [5]). 

The crucial point in studying quantum correlations in the system of unstable 
particles, say K^K^, is the choice of a model describing time evolution of the 
system under consideration. It is obvious that this choice depends on what 
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physical aspects we want to neglect and what we want to focus on. For exam- 
ple, if we want to focus our attention on decay of particles governed by the 
Geiger-Nutall law neglecting the evolution of decay products and other physi- 
cal processes hke e.g. decoherence, we can use the standard Weisskopf-Wigner 
approach [20,21]. The price we must pay for this choice are some ambiguities, 
due to the non-Hermitian Hamiltonian of Weisskopf-Wigner model, when we 
want to consider a series of consecutive measurements. However, if we mainly 
want to take into accout decoherence we can do our calculations in the frame- 
work of many models considered in the literature (see e.g. [22,23,24,25]). The 
common feature of all these models is taking into account only the two-particle 
sector of density matrix describing the state of the K^K^ system, neglecting 
the one-particle and zero-particle sectors which arise during the time evolution 
due to the decay process. 

Recently, we have introduced a model [26] of time evolution of a neutral kaon 
in its rest frame, which allows one to describe the K^K^ system in more realis- 
tic way. In the framework of our model, based on the theory of open quantum 
systems (see e.g. [27,28]), we can take into account that the state of the sys- 
tem evolves from a two-particle state to the zero-particle one through states 
being a mixture of all possible states. This model is based on the following 
assumptions, (i) A neutral kaon is treated as an open quantum system, and 
the decay products are regarded as a part of the environment, (ii) The system 
can be found either in a particle state, or in the state of absence of the particle, 
which we call (in analogy to the quantum field theory) the vacuum state, and 
denote by 1 0) . Therefore the Hilbert space of the neutral kaon-vacuum system 
is a direct sum of one-particle Hilbert space and one-dimensional space Tig 
spanned by |0). (iii) The superselection rule prohibiting the superposition of 
the particle and vacuum holds, (iv) The time evolution of the system, consis- 
tent with the Geiger-Nutall law and allowing for CP symmetry violation, is 
given by a family of completely positive trace preserving maps forming a one- 
parameter dynamical semigroup. Complete positivity implies that the time 
evolution of a state of the system can be written in the operator-sum (Kraus) 
representation [29] (which immediately allows one to find the time evolution 
of non-interacting particles) 

i 

and the trace preservation requirement leads to the condition 

Y,Ei{t)m-^- (2) 

i 

As a remarkable side effect of complete positivity of time evolution of the sys- 
tem we obtained the upper bound for the decoherence parameter (see [26]). 
Here we would like to point out that we used the dynamical semigroup ap- 
proach for the entire evolution of the system, not only for the description of 
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its decoherence, as was done in [22,23,24,25]. 

Tlic paper is organized as follows. In Section 2 we generalize our description 
of time evolution of an unstable particle in its rest frame to the case of an 
arbitrary reference frame. Then, in Section 3, we show how the evolution of 
the system of two unstable non-interacting particles can be found when we 
know the Kraus operators governing the evolution of its components, both 
for distinguishable and identical particles case. Finally, in Section 4 we apply 
these results in calculation of some quantum correlation functions for K^K^ 
system in the singlet state. And we show that the obtained results are exactly 
the same for kaons treated as distinguishable or identical particles. 



2 Time evolution of an unstable particle in an arbitrary reference 
frame 



Let us assume that we know the Kraus representation of time evolution of an 
unstable particle in its rest frame in which the four- momentum is A; = (m, 0). 
Let us answer the question, how this evolution looks like from the point of 
view of an observer in the laboratory frame, in which the particle has four- 
momentum k = {k^, k). Let Lk be the Lorentz boost from the rest frame to the 
laboratory frame, i.e. k = Lkk. When we denote the inner degrees of freedom 
(e.g. strangeness, bottom) as s, than the state of particle with four-momentum 
k can be denoted as \s,k) (in this notation, |0) = |0, 0)). Since we deal only 
with (pseudo) scalar particles (e.g. kaons), \s, k) = U{Lk)\s, k), which implies 
that the density operator in laboratory frame has the following form 

t'^it) = U{Lk)g{T')U\Lk) 

i 

where the proper time in the rest frame of the particle is denoted by r'^ because 
it is convenient to express it by means of k and the laboratory time t: 

= ^ = — (4) 



and ^'^^(r*^) = U{Lk)Ei{r^)W (Lk). By means of the decomposition Ei{T^) — 
Y.s,s' El^' k){s', k\ (see Appendix A), we get 

E^'^\T') = j:Er'{r%s,k){s',k\. (5) 

s,s' 

It means that the matrix elements of E^^'^^(^t^) in the basis {\s, k)} are exactly 
the same as the matrix elements of Ei{T^) in the basis {\s,k)}. Note that this 
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holds only for (pseudo) scalar particles. 

Taking into account the form of the Kraus operators given in Appendix A, 
the evolution (3) can be easily extended to the case in which the particle state 
is the superposition (or mixture) of a few different momentum eigenstates. In 
this case the space of states of the system is spanned by the orthogonal vectors 
\s, k) and vacuum, but now the four-momentum k can take an arbitrary but 
finite number of different values. We denote the set of admissible four-momenta 
by Q (for example, in the case of two identical particles discussed in the next 
section, Q consists of only two elements Q = {p,q})- The time evolution of 
the density operator describing such a state has the following form 

5 

k&Q i=l 

where the Kraus operators Eq = E^^ = |0, 0)(0, 0| and E^^^^It'^), are given in 
Appendix A. Note that when Q consists of only one element Eq. (6) reduces 
to (3). 



3 Two unstable peirticles as open quantum system 

Let us consider two unstable particles A and B (this approach can be eas- 
ily extended to the multiparticle case), with four-momenta p and q in the 
laboratory frame, respectively. The space of states of the whole system is 

n^iUA® no) (Hb e Ho) 

= (Ha ® Hb) e (Ha ^Ho^Ho® Hb) © {Ho ® Ho), (7) 

so one can easily see that the system can be either in two-particle state 
{Ha <8) Hb): or in one-particle state {Ha ®Hq®Hq® Hb), or in the vac- 
uum state {Hq ® Ho)- Such a system is an open one, so its evolution is repre- 
sented by a completely positive map and can be written in the operator-sum 
representation 

Q'X%{t)- E E^rW^l^(0)(^r(^))^ (8) 
keQ,k'eQ' i,j 

where, under a reasonable assumption that the particles do not interact, we 
define 

E^f(t) = if "^(r'^) ® Ep^^{T^'). (9) 

E^^^^{t) are the Kraus operators governing the time evolution of a one-particle 
system, and i, j = 0, . . . ,5 (see Appendix A), and Eq^^^{t) = Eq. The normal- 
ization condition J2keQ,k'eQ'T.i,j{Eff {'t)VE^j'' {t) = 1<8)1 is fulfilled, which 
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follows from corresponding normalization conditions for one-particle Kraus 
operators (A. 2). It is clear that for distinguishable particles we have Q = {p} 
and Q' = {q}. However in the case of identical particles one does not know 
which of them carries which four-momentum or flavour, thus the space of the 
whole system must be symmetrized, and Q — Q' — {p, q}. Let us introduce 
the permutation operator V such that V{\s, k) (8) k')) — k') (8) |s, k); then 

Vg'XB{t)V = g'X^'sit)- (10) 
Moreover any observable O must also preserve symmetrization, i.e.: 

vor^o. (11) 

One can easily verify that if ^^^(0) fulfils the condition (10), than also Q'^it) 
obtained from (8) fulfils it, so the time evolution of identical particles is gov- 
erned by the same evolution law (8) as the time evolution of the distinguishable 
ones. 

In EPR-typc experiments it is of great importance that the measurements 
performed by distant observers, say Alice and Bob, must be local. The local- 
ity means that the observables are restricted to some specific region usually 
interpreted as the region of detector (see e.g. [30]). It can be achieved exper- 
imentally by assuming, for example, that Alice's detector can register only 
particles with four-momentum p and Bob's only those with four-momentum 
q (directions of momenta p and q must be sufficiently different) ^ . Thus, only 
the particle carrying four-momentum p and q, respectively can reach Alice's 
and Bob's detectors. 

The space of states of a two-particle system is the tensor product (7Y^ © Tio) ® 
{Hb © Ho) (see (7)). In the case of distinguishable particles T-Lj, is spanned by 
and T-Lb by so Ahce's detector registers only the particles 

with states from subspace Tiji^ © TIq and Bob's only those from Tis © 'Hq. On 
the other hand, in the case of indistinguishable particles Ha and Hb must 
be identical and both spanned by the same set of vectors {|s,p), The 
physical Hilbert space is symmetric subspace of H, i.e. |(l-|-P)7i, therefore 
one cannot associate specific one-particle Hilbert space with a given observer. 



4 Quantum correlations in the neutral kaon system 

Let us consider two neutral kaons in a given initial state Q^^{^)-, and two 
distant observers, Alice and Bob, in the same laboratory frame. Alice can 

^ Actually the sharp momentum states cannot be achieved because of the uncer- 
tainty principle and finite volume of detector. 
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measure the flavour of kaons only with four-momentum p and Bob only those 
with four-momentum g, therefore their observables commute, i.e.: 

[>,S«] = 0. (12) 

Hereafter, we omit the superscript 'lab', because in the sequel we consider 
only density operators as seen from the laboratory frame. 



^.1 Distinguishable case 



In this subsection we treat the kaons as distinguishable particles, as it is usually 
done to simplify calculations. The Hilbcrt space of the first kaon is spanned 
by set of orthonormal vectors \K^,p)^ |0)0)} ^"^^ the Hilbert space 

of the second one by q), \K^-, q), |0, 0)|. 

Suppose that at time tA Ahce measures an observable — A<Si 1 with the 
spectral decomposition J2a<^(^a® 1)- The density operator just before the 
measurement is 

QAsitA) =T.E!!{tA)QABmmtA))\ (13) 

(cf. (8)), and when the outcome of the measurement is a, the state reduces to 

QAB\tA) = —J—. , (14) 

Pa[tA) 

where Pa{tA) is the probability of measuring a at time Ia- Next, at time ts, 
Bob performs the measurement oi B'^ — 1 <S>B with spectral decomposition 
X^5 6(l(8)nb). Just before his measurement the state is 

QAsitB) =Y.EUitB -tA)gAB{tAm'{tn -tA)y. (15) 
ij 

The conditional probability that Bob's outcome is b provided that Alice's was 
a is 

p^„itB) = TV {(1 ®n,) ^^(^b) (1 m,)} . (16) 
By means of (13)-(16), the joint probability 

Pab(tA,tB) = Pa(tA)Pb\a(tB), (17) 

that Alice's and Bob's outcomes are a and b, respectively is given by the 
formula 

PabitA, ts) = TV{^^(0) [£{Ef (r^))tn„^f (r^) ® ^{E] {t^b^.E] (r|)] }. 

i j 

(IS) 
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Then the correlation function between the outcomes 



CAvBi{tA,tB) = Y^ab Pab{tA,tB), (19) 
ab 

takes the form 

CArB.{tA,tB) = Tr{^^B(0) [A{rl) S(t|)]} , (20) 
where yl(r^) = E.(£^f (r^))^^^f (r^), and 5(rl) = E.iE] {t%))^BE^, {t%) . 

Now we calculate explicitly a few correlation functions and probabilities in 
the system of neutral kaons in the pure entangled state J^'~' —1 (produced 
through the reaction e+e" 0(1020) K^K^) 

1^) = i= [\K\p) ® \K\q) - \K\p) ® \K\q)) (21) 

to show how our model works. We do not neglect CP violation and decoher- 
ence, despite the fact that CPLEAR experiment at CERN (where K^K^ pairs 
are produced in the pp collider) is not sensitive to CP violating effects. Let us 
begin with the strangeness correlations. The strangeness operator takes the 
form 

S^= \K'^,k){K\k\ - \K\k){K'^,kl (22) 

where k = p,q, and — 1, B'^ = 1 (gjS'^. It is easy to show that the 
corresponding correlation function has the following form 

CsMtA,tB) = -^-i^[e-(^+^)(-A+-^)cos(AmAT) 

- l^l (e-^«^^-^^^^ + e-r^-^-^^-s)] , (23) 

where and Tg stand for the proper times, Ar = t^— r^, Sl — 23?(e) / (l+|ep) 
(e is a small complex CP-violation parameter), T = ^(Ts + Ti,) (Ts and F^, are 
the decay widths of short and long living states of neutral kaon, respectively). 
Am = rriL — ms {rns and ttil are masses of short and long living states of 
neutral kaon, respectively) and A is a decoherence parameter, representing 
interaction between one-particle system and the environment. On the other 
hand, the strangeness operators defined above have three different eigenvalues 
±1 and 0, but the observables considered in Bell-CHSH inequalities have only 
two different eigenvalues ±1. Such a dichotomic observable is, for example, an 
observable answering the question whether one registers a kaon (anti-kaon) 
(then the result of the measurement is -|-1), or not (then the result is — 1). 
Such a case was investigated in [3]. Denoting this observable as {D-), in 
one-particle case we have 

= ±\K^,k){K'^,k\ T \K°,k){K^,k\ - |0,0)(0,0|. (24) 
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It is easy to show that in the case when — D\ ® 1 and W — 1 (8)D+, we 
get 



1 + 5l r _(r+A)(.^+.|) cos(AmAr) 



2(1 -5l)^ ' 
+ f e-(^+^)^^ cos( Amr^) + e-^^+^^^B cos(Amr|)) , (25) 

and when k' ^D\®\ and S« = 1 ®Dl, we have 

Cd^di (^a, ^b) = 1 + e-(r+^)(^A+^^) cos(AmAT) 

_ ^^Y^yj (e--^r« + e--Ar^ _25^ e-^^Cr+A) cos(AmT5)) 

- , V , fe-^^rs + e-^B^^ +25^ e-^s(r+A) cos(AmTl)) 
2(1 + ^ 

+ 1 (^e--Ars-r^r^ +e-"^r^-^s^s) . (26) 



Of course, we could find the above quantum correlation functions directly 
from the definition (19) finding appropriate probabilities. For example, one 
can easily find the stangeness correlation function Cspsi{tA,tB) knowing that 
(i) the probabihty that Alice's detector registers at tA and Bob's K'^ at 
ts is 



PK0,K0[tA,tB) = 7- 

o L — Ol 



- 2 e-(^+^)(^A+^l) cos(AmAT)l , (27a) 



(ii) the probabihty that Alice's detector registers at tA and Bob's K'^ at 
ts is 



PKO,Ko{'tA, ts) — o 



SI + Sl 



- 2 e-(r+^)(^A+^^) cos(AmAT)l , (27b) 



(iii) the probabilities that Alice's detector registers at and Bob's at 
ts, and that Alice's detector registers at tA and Bob's K'^ at are 

PKO,Ko{tA,tB) ^ PKO,Ko{tA,tB) 

^ 1 [e-rsr^-r^r^ ^e-^^^s-r^^^ +2 e-(r+^)(^A+^^) cos(AmAT)" . 

(27c) 
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In [3] the correlation function C^jp^)? (tyi, t^) and probabilities PKO,Ko(tA,tB), 
PKO,Ko{tA,tB) and PKO,Ko{tA,tB) wcrc found under condition 5^ = A = 0, 
i.e. without CP-violation and decoherence. Of course, when we put (5^ = A = 
in (25) and (27) we arrive at Bertlmann's results, if we only take into account 
differences in conventions. 



4-2 Indistinguishable case 



Now let us consider the same situation as in the previous subsection but with 
more realistic assumption that kaons are indistinguishable particles. Suppose 
that Alice measures at time and Bob measures 13'^ at time ts- Both of 
the observables must fulfil (11). After analogous calculations as in the previous 
subsection, using (12), we get 

CA.B.itA,tB) = Tr [^^(0)i^(r^)^^(rl)] , (28) 

where 

Mr'A) = Y:imtA)yA^EfJ{tA) and B'^{t],) = Y.iE^jitB))'B'^E^!{tB) ■ 

ij ij 



Now let us calculate the same correlation functions as before. First, we have to 
note that in the case of indistinguishable particles the initial state has different 
form 

IV^) = ^ {\K',p) \K', q) + \K', q) ® \K\p) 

- \K\p) ® \K\q) - \K\q) ® \K\p)) . (29) 

Second, we have to construct observables that answer the same questions as 
the observables used in the case of distinguishable kaons. The strangeness 
operators take the form 

= s''®i+i (30) 

where k takes the value p or g, and was defined in the previous subsection. 
Observables D\ cannot be constructed in analogy to (30) by means of (24) 
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and symmetrization ^ . Now, the form of is 

Dl^2 k){K^, A;| 1 + 1 k){K^, k\) 

- 1^1 -\K\k){K'^,k\ (g) \K^,k){K^,k\. (31) 

It yields +2 when both of the particles are kaons with four-momentum k, +1 
when one of the particles is a kaon with four- momentum k, and —1, when there 
is no such kaon. It is not dichotomic, but it is not a problem because in state 
(29) the probabiUty of measuring two kaons with the same four-momentum 
equals zero. D'^ takes the analogous form to l)^ 

Dt^2(\K^,k){K^,k\®l + l(^\K^,k){K^,k\) 

-l0l-\K^,k){K^,k\ ® \K^,k){K^,k\. (32) 

It is easy to check that the correlation functions Cgpgq{tA,tB), Cf^v fyq (^Aj^b) 
and CfjPjjq {tA-ts) arc exactly the same as in the distinguishable particles case 
(23), (25) and (26). It is quite remarkable that for the observables considered 
above, it does not matter whether we treat kaons as indistinguishable particles 
or not, at least in the singlet state. Of course, it is easier to carry out all 
calculations on the assumption that kaons are distinguishable particles, as it 
is usually done. 



5 Conclusions 

We have shown that having the Kraus representation of time evolution of 
an unstable particle in its rest frame it is possible to find the evolution of 
the particle in an arbitrary reference frame. Moreover, we have also shown 
that taking into account the form of Kraus operators one can extend the 
time evolution of the system to the case in which the state of the particle is 
the superposition or mixture of a few different momentum eigenstates. Next, 
we have found the time evolution of two non-interacting unstable particles, 
distinguishable as well as identical ones. Finally, we have applied these results 
in calculation of some quantum correlation functions for K^K^ system in 
the singlet state assuming CP- violation and decoherence. And it turned out 
that the results are exactly the same either we treat kaons as distinguishable 
particles or identical ones. Therefore, one can neglect the fact that kaons are 
identical particles and treat them as distinguishable ones, at least in the cases 

^ When we calculate the spectral decomposition of, say, ® 1 + 1®!)^, where 
= '2,\K^ , k){K^, k\ — 1, we find out that it has eigenvalues equal ±2 and 0, so it 
does not answer the question, whether the particle is kaon carrying momentum k, 
or not. 
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considered in this paper. It is still an open question whether the statistics of 
the particles plays no role in the general case or not. 

We would like to point out that all results presented in this paper will be 
valid also for B-mesons after appropriate change of notation, because kaons 
and B-mesons evolve according to the same scheme. 
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A Time evolution of 



Let us assume that the space of states of the system is spanned by the or- 
thogonal vectors: vacuum |0, 0) and \s,k), where s denotes inner degrees of 
freedom and four-momentum k belongs to the finite set Q of admissible four- 
momenta. The Kraus operators Ef{t) = J2s,s' E-'"^ {T'')\s,k){s',k\, describing 
evolution (6) of the kaon which state can be a superposition (or mixture) of 
different momentum eigenstates, have the following form 



E',{t')^Eo=\0,0){0,0\, 



_^ g-rHA+2imi+rz,)/2^ 



X [\K\k){K\k\ + \K\k){K\k\) 

_^ 1 '^''('*^+2«"^S+rs)/2 _ g-Tfc(A+2imi+ri)/2^ 

X (\±l\K\k){K\k\ + ^-^\K\k){K\k\ 



1-e' 





e2 


- / 







l_e-'=r, _^2 



1 + e' 

_ g-T'=(r+A-iAm) 



v 



1 _ e-^'^rz 



X 



^ ■0,0){K^,k\ + -^\0,0){K^,k\) , 



1 + e 



1-e' 



(A.la) 



(A.lb) 



(A.lc) 
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1 + 



'1 - e-^'r^ _ g-r'^CA-iAm+r) 



^ 2(1 - e-^'^r^) V 1 + e 

1 - e-^'r^ _r ^ e-r'=(A-iAm+r) x 



1-e 



|0,0)(i?°,^| 



0,0)(i\:°,A;| 
(A.ld) 



1 + e 



1-e,.-. 



4fc(r'=) = ]-e-^''^^l^^l^^^^(\K\k){K^,k\ + A;| 



1 



1 + e' 



(A.le) 



(A.lf) 



where = t/j^. e is a small complex CP- violation parameter, = 2K(e) /(1+ 
|ep), and are the decay widths of K*^ and (short and long living 
states of neutral kaon), respectively, F = |(F5 + F£,), ms and rriL are masses of 
Kg and K^, respectively. Am = — m^, and A is a decoherence parameter, 
representing interaction between one-particle system and the environment. In 
comparison to [26] we use different, more convenient set of Kraus operators, 
leading to the same evolution. It is easy to check that the normalization con- 
dition 

5 

eIEo + ^ ^(if ^"'^(r'^))tif i^>^(r'=) = 1 (A.2) 

fcGQ «=1 

holds. In the case of indistinguishable particles the Kraus operators are exactly 
the same, but we additionally must sum over all admissible k. 
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